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ABSTRACT 

A method of analysis capable of predicting accurately the fracture 
behavior of a unidirectional composite laminate containing symmetrically 
placed buffer strips is presented. As an example, for a damaged graphite/ 
epoxy laminate, the results demonstrate the manner in which to select the 
most efficient combination of buffer strip properties necessary to inhibit 
crack growth. Ultimate failure of the laminate after crack arrest can 
occur under increasing load either by continued crack extensiov? through 
the buffer strips or the crack can jump the buffer strips. For some 
typical hybrid materials it is found that a buffer strip spacing to width 
ratio of about four to one is the most efficient. 

INTRODUCTION 

One of the major difficulties ir. designing an advanced composite 
structure such as an aircraft to comply with current safety regulations, 
is meeting the damage-tolerant (fail-safe) requirements. One very 

1 This work was supported by the Fatigue and Fracture Branch, Materials 
Division, NASA-Langley Research Center unc.er Grant NSG-1297. 




promising method of constructing s damage-tolereut composite laminate Is 
to use hybrid, embedded stringers (buffer strips) es e creek errest 
mechanism. A typical laminate Is shown In Figure 1 and the geometry 
assumed for the present study Is given In Figure 2. Two fundamental dif- 
ferences are seen between the real construction and the model; first the 
model is assumed to consist of only xero degree (parallel to the load) 
fibers and second, it contains an initial central crack between two 
buffer strips and two half-planes rather than a periodic array of buffer 
strips. It is felt that much of the characteristic behavior can be re- 
presented by the unidirectional laminate, as a dominant portion of the 
load is carried by these fibers. A primary function of the angle plies 
in Figure 1 is to prevent longitudinal matrix splitting in a brittle 
matrix such as epoxy. This is accounted for to some degree in the present 
solution by allowing the matrix to support large strains without splitting. 

This work is an extension of the studies presented by the authors in 
references 1, 2 and 3 and is the latest solution developed in an attempt 
to understand the damage tolerant behavior of a buffer strip laminate. 

The intent is to be able to estimate the remote stress required to fail 
the hybrid unidirectional laminate of Figure 2. The fibers and matrix 
are assumed to be linearly elastic and the failure criterion is simple 
tension failure of the fibers. The classical shear-lag model is used to 
represent the shear stress distribution between adjacent fibers. From 
the previous work^ it is known that for a single material laminate, with- 
out matrix yielding and splitting, the most highly stressed fiber is the 
first unbroken one directly in front of the notch. The significant 
question in this study is, if the first fiber in front of the notch breaks 
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at a given applied strees will che next fiber require a higher or lower 
applied stress co also fail or will a fiber break ac some ocher locadon 
at an even lower screes? Thac is, is che crack growch stable or unstable, 
and how does this behavior depend on tsaterials and geonetry? The shear- lag 
model and che assumptions and simplifications made in che analysis are 
somewhat restrictive but Che ability of chese simple models to represent 
actual laminate response has been found* to be very good, it is then felt 
chat che results given in this paper are a good indication of the behavior 
of a buffer strip laminate. 

The initial studies using the shear-lag model to analyse notched uni- 
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directional laminates were given by Hedgepeth and Hedgepeth and Van Dyke * . 

The work of 1, 2 and 3 extends these methods up to the present treatment. 

Experimental investigations concerning buffer strip laminates are dis- 
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cussed by Eisenmann and Kaminski ; Hess, Huang and Rubin ; Vvery and 
Porter^; Verette and Labor*^; and Poe and Kennedy**. Because of the 
limited space allowed for this presentation much of the background de- 
tails and development must be referred to these papers. 

FORMULATION 

The fundamental solution needed in the analysis of this problem, is 
the case of a unidirectional half-plane with broken fibers and matrix 
splitting as shown in Figure 3. This basic solution will be developed 
first, and then by taking appropriate combinations of particular forms of 
this result, the complete solution will be presented. 

A unidirectional array of parallel fibers with an arbitrary number of 
broken fibers in the form of a notch and a longitudinal split in the matrix 
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is shown In Figure 3. The laminate is subjected to prescribed sheer 

stresses T^(y) along the free edge end T fe (y) along the split, end a remote 

uniform tensile strain in the axial direction. Fiber breaks occur along 

the x-axis (axis of symmetry) and, since the loading is e y e trie, only 

the upper half of the laminate is considered in the analysis. 

The fibers are taken to be of much higher strength and extenslonal 

stiffness than the matrix and all of the axial load is assumed to be 

carried by the fibers with the matrix transferring load by shear stresses 

4 

as given by the classical shear-lag assumption . The axial fiber stress, 
o n (y) , and matrix shear stress, T Q (y) , are then given by the simple re- 
lations 

dv (y) Gw 

°n (y) " E F~d7~ * “ d T n (y) " T [v n (y) ' V n-l (y)1 ’ (1) 

where, v (y) is the axial displacement of the fiber n at the location y, 

E F is the Young's modulus of the fiber, is the equivalent matrix shear 
modulus and h is a shear transfer distance. Because of the interference 
between fibers it is unlikely that will be the homogeneous matrix shear 
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modulus or h the actual fiber spacing. It is pointed out by Goree that 

these values can be determined experimentally for a given laminate. 

12 

Batdorf also discusses this question in considerable detail* 

By virtue of the shear-lag assumption the longitudinal and transverse 
equilibrium equations become decoupled and the fiber axial displacements 
and stresses can be obtained without solving the transverse equilibrium 
equation. Therefore, only the equilibrium equation in the longitudinal 
(axial) direction will be considered. With reference to the free-body 
diagram of a typical fiber-matrix region shown in Figure 3, the equilibrium 
equations in the longitudinal direction are given by 
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t dy * 1 
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x (y) - * m (y) " o » for fiber 0, 

for fiber n. 


K do (y) 

r-$r+ w*>-\w • 


d °nw ( y> 

t dy 


+ x K (y) ~ ^(y) ■ 0 , for fibor HW when y <, * » and 


d °NW-l (y) 


t jy + T NW-l (y) ** T b (y) " 0 » for £lb * r NW * 1 * l#0 y - 4 ’ (2) 

Using the s tress -displacement relations , Equation (1) , in the above equilib- 
rium equations, the following set of differential-difference equations is 
obtained: 




V 

+ V 

dy 2 

w 

d^v 

!L + v 

G M t 

. 2 + v i 

dy 

Vr h 

h2 v 

a V NW+1 

V 

dy 2 

w 

h2 v 

NW+1 

V 

dy 2 


+ v. - v - x (y) , 


0 , 


- v, 


NW 
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Noting the coefficient of the second derivative term in the above 
equations, the following changes in the variables are suggested: 


let y ■ 


w 




d y n 

a a * E_ -3 — , and 
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Algebraic manipulation then gives 
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^n 

c n " — * T n 


Vr 

^ (V n-Vl )ttd 4 “^| 


w 

V 


8 , 


(5) 


where n . 8 , o end V (n) are non-dimensional . 
n n 
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By making use of Fourier transform techniques the resulting dif- 
ferential-difference equilibrium equations can be written in the form of 
a single differential equation given by 


i ^y(n» 9 ) . ^(n.e) * t a (n) cos(y) + <n-B>[g(n) - t b (n)]F 1 2 3 , (6) 

dri 

where, 

V n (n) " f ^ ^ (n,6) C08[(n + I )e]de • " 2[1 - cos(6) ] * 4 sin 2 (|) , 


V n) " 


EFth t (y) t b (y) 

-• V"> ’Na^ — 


»- 

F 2 - cos [ (NW + i)6] - cos [ (NW + |)6] , 


< t \ - g> 


1 for n <. 8 
0 for n > B , and 


" V NW+1 “ V NW 


The solution to the problem of vanishing stresses and displacements 
at infinity and uniform compression on the ends of the broken fibers will 
now be sought. The complete solution is obtained by adding the results 
corresponding to uniform axial strain and no broken fibers to this 
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solution. As before the solution to Equation (6) satisfying vanishing 
stresses and displacements reduces to solving a set of linear algebraic 
equations, in terms of the unknown Fourier constants , given by 
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+ ■- / cos (4) cos [(n + 4)0]/ .“ 6t x (t)dt d0 

* o 1 * o * 

+ — / F 2 cos [ (n + 4)0]/ s“ fi# (g(s) - t. (s) } ds d0 » 1, (7) 

* o 2 o ® 

for all broken fibers, l.e., n ■ 

The displacement of any fiber n at n is then given by 

V (n) - “ /* e“ 6n I B cos [ (N* + m + 4)0] cos!(n + i)0]d0 

n w o m-i m 2 1 

- :/' SS*M2 1 f m D(«,n,t) t (t)dt cos[(n + ^OjdO 


* o 6 o 


» V 2 0 


- i ! D(6,n,t){g(t)-x b (t)>dt cos [(n + |)0]d0 , (8) 


where, D(5,n,t) * e 


o 


. i „ ,.i 

-0 |h“t| 


- e 


-6 (n+t) 


SYMMETRIC BUFFER STRIP LAMINATE 

Since the laminate shown in Figure 2 is symmetric about the x and 
y-axes , only the upper right quadrant will be considered. Figure 4 shows 
the three distinct regions of the laminate. Regions I and II are finite 
width unidirectional strips with broken fibers subjected to remote tensile 

stresses o 1 and o^and varying shear stresses along the free edges. The 

00 00 

solution of these two regions can be obtained by setting the split length 

equal to infinity in the basic solution obtained in the previous section. 

The region III is a unidirectional half-plane subjected to uniform remote 

tensile stress o 111 and varying shear stress t^^Cy) along the free edge, 

" b 

the solution of which is obtained by setting the split length equal to zero 
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In the untie solution of ths previous ssetion. Thus tha solutions for all 

ths thrss regions are known for given applied shear end axial straasas. 

Where these regions ere joined together the sheer stress is unknown. 

But from equilibrium, the sheer stresses on each of the adjacent regions 

must be equal at their respective interfaces. Further, as the sheer 

stress is directly related to the distortion of the matrix from the sheer- 

leg assumption, it follows that these stresses must be proportional to 

the difference in the displacement of the adjoining fibers of the adjacent 

regions . These conditions, along with the stress boundary conditions on 

the broken fibers in regions I and II, are used in obtaining the solution 

for the entire buffer strip laminate. The superscripts I and II indicate 

/ i2 

the variables in region I and II. Further, sod ate the 

ratios of and h for interfaces I and II, respectively. Denoting 

f 1 (t ) - T I (n) - g I (n) , f n U) - "“(5) - g U (0 , 
a b 

F 1 » cos [ (NW1 + y)6] - cos [ (NW1 + -|)e] , 

F 11 - cos [ (NW2 + y)0] - cos [ (NW2 + |)0] , and 
C(k) - cos[(k + ^)9] 


the governing equations for the buffer strip laminate can be given as 
follows : 


t 

71 o 


2 * 


( M. 

I B 1 C(m)6 - F 1 / e‘ 6t f X (t)dt 
m*0 m o 


C(n)d6 - 1 , 


( M 


I B 11 C(N* + m)5 + G., / e 1 T X (t)dt 

m 2 lr 


tn«0 


R o 
1 


_II -os ,11, x. 

- F ; e f (s)ds 


C(j)d6 - 1 , 


(9) 
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(12) 
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12 R 2 5 o 1 * 
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F II C(NW2)f II (s) + C 2 (0 )tJ X (s) } ds ) d6 


Mr 
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T XI (0 • / ff | - 2 l B XI C(N* + n)C(NW2)e‘ 65 


m«0 
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+ Gi2 ms r P ( i ,t/R 1 .o7 2 ct)dt 

- r^F 1I C(NH2)D(6,E,s)f II (s) + 0 23 ^-^-D(e 2 ,5,s)T 21 (s)j (tajae. 


where, and are the number of broken fibers in 1 and II, 
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- fi/l^ » « 2 - /Rj , G u • (G^/h) 11 (WGj,) 1 » 

C i2 - oyh) 12 (h/G m ) IX Md - (G^/h) 11 (WGg) 111 . 

The above governing equationa are of the same fora aa thoae obtained In 

3 

the caae of a • ingle buffer atrip laminate Dharani end Goree , except 
that this problem has an additional Integral equation due to the finite* 
ness of the center panel. 

SOLUTION 

The above six equations (9-15) contain the unknown Fourier conatants 

B X and B XX and the unknown functions g X (q), g XX (0» T X (h)» end T XX (0« 
m m • o 

The solution is developed by representing the integrals constaining the 

3 

unknown functions using a Gauss-Laguerre quadrature formula and reducing 
the six equations to a single system of equations having as unknowns the 
Fourier constants and the values of the unknown functions at specific points 
(quadrature points) . 

For any continuous* integrable function the Gauss-Laguerre quadrature 
formula gives 

eo -»X 

f f(x)dxsl w.e 1 f(x.) (15) 

o i-1 

where is the i th zero of the Laguerre polynomial, L^x^ , and w^ is the 
corresponding weight function given by 

- x 1 /[(K+l)L K+1 (x i )] 2 . (16) 

For the results presented in this paper, forty-five terms (K“45) were 
taken to represent each of the four unknown functions. Computation time on 
the Clemson University IBM 3081-K computer was about two minutes for a 
typical geometry. 
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Figure 5 presents results corresponding to Initial crack growth in 
region 1, crack arrest at tha Interface, creek growth In the buffer strip 
and subsequent laalnate failure. In these results all fibers are of the 
sane cross-sectional area and In all cases the buffer strips are ten fibers 
wide end are thirty fibers apart. Two buffer strip naterlals are con- 
sidered, each with the sane modulus but with different ultiaate stresses as 
shown in Figure 5. Material 2 has properties close to that of 8-glass 
end the parent laminate Is graphite/epoxy. The solid line in Figure 5 
represents the renote stress required to Initiate crack extension, (fall 
the first unbroken fiber in front of the notch, fiber A). The renote 
stress required to fall the laalnate catastrophically, (fail the first 
fiber it plane III, fiber B) is given by the broken line In Figure 5. Both 
these stresses are functions of the initial crack length end decrease with 
increasing length. Results for an all graphite /epoxy laminate are also 
given. The crack growth takes place by breaking consecutive fibers fron 
the crack tip to the Interface. Then, depending on the stress level 
required to run the crack to the Interface and depending on the buffer 
strip material, the crack may arrest. Both buffer strip naterlals require 
an Increasing stress to continue the crack growth in the buffer strip, 
although material 1 will arrest a crack only if it initiates under fairly 
low load, l.e. initially close to the interface. For the particular lamina 
of Figure 5, all fibers in the material 1 buffer strip fall before fiber 
B attains its failure stress, whereas for material 2 fiber B falls when 
there are still some fibers left unbroken, l.e., the crack jumps the buffer 
strip. 


In Figure 6 the effect of buffer strip width on creek growth for e 
fixed epeclng between buffer stripe of thirty fibers Is given. The 
ultlaete f el lure stress of the lealnete es e function of buffer strip width 
Is plotted In Figure 7. From Figure 7 It Is seen thet for aeterlel 1 the 
optlaua buffer strip width Is shout 3*4 fibers end for aeterlel 2 , ebout 

8 fibers. Additlonel results Indicate thet one nay think of Individual 
fibers as groups of fibers end Figure 7 then lsplles thet, for e (graphite, 
S-glass) /epoxy hybrid lealnete, the optlaua espect ratio should be about 
four to one. 
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Legends for Illustrations 


Figure 1. 
Figure 2. 
Figure 3. 
Figure 4. 
Figure 5. 
Figure 6. 
Figure 7. 


A Typical Buffer Strip Laminate. 

Geometry of s Sysmstric Buffer Strip. 
Unidirectional Half-Plane with Broken Fibers. 
Three Regions of the Buffer Strip Laminate. 
Effect of Buffer Strip Width on Crack Growth. 
Failure Stress as a Function of Crack Length. 
Ultimate Failure Stress vs. Buffer Strip Width. 
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